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Abstract 

We reexamine several issues related to the physics of scaling in electron scattering from nuclei. 
A basic model is presented in which an assumed form for the momentum distribution having both 
long- and short-range contributions is incorporated in the single-particle Green function. From this 
one can obtain saturation of nuclear matter for an NN interaction with medium-range attraction 
and short-range repulsion, and can obtain the density-density polarization propagator and hence 
the electromagnetic response and scaling function. For the latter, the shape of the scaling function 
and how it approaches scaling as a function of momentum transfer are both explored. 

PACS numbers: 24.10.Cn; 25.30.-c 
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I. INTRODUCTION 



Scaling phenomena as realized in electroweak interactions with nuclei at intermediate-to- 
Ugh energy, have .eve.al facet, include .caUng of the a.t kind Q (independence of tlte 

momentum transfer q) and second kind [2|, |3[ (independence of nuclear species, for instance, 
as characterized by the Fermi momentum kp)- In addition, universality for different reaction 

nn 

channels (longitudinal, transverse, etc.) has been called scaling of the zeroth kind [3|, 
while universality in the isoscalar and isovector channels has been called scaling of the third 
kind [5]. All four types of scaling are reasonably well respected by data at sufficiently high 
energies, namely away from threshold and for momentum transfers typically twice the Fermi 
momentum or larger, although there are observed to be scaling violations and their origins 
provide interesting insights into the dynamics of the scattering processes. For example, 
the transverse EM response is known to involve both the familiar one-body currents and 
two-body Meson- Exchange Currents (MEC) which do not scale in the same way {g-S. 10. 



In the present study we focus on one particular aspect of scaling, namely, scaling of the 
first kind. Our motivation is to explore the interconnections between the strong interaction 
dynamics of a representative NN potential that is chosen to provide the correct binding 
energy and saturation density of nuclear matter, on the one hand, and a corresponding Green 
function that is made to be consistent with those properties of nuclear matter, on the other. 
Having such a Green function one can immediately obtain the density- density polarization 
propagator and hence the longitudinal response Rl and scaling function Fl. For brevity 
these are called simply R and F in the rest of this work. In contrast to the usual approach, 
in the present study we "work backwards" assuming a form for the momentum distribution 
n{k) and, given this, obtaining the corresponding energy per particle as a consequence. In 
particular, the chosen momentum distribution is taken to have both long-range contributions 
(those below and slightly above the Fermi surface) and short-range contributions which give 
rise to a tail that extends to high momentum. 

To be able to carry out this study with much of the development still analytic we restrict 
our attention to the non-relativistic situation and assume a translationally invariant (infinite, 
homogeneous) many-body system of point nucleons. Our goal is not to provide a detailed 
numerical study of scaling phenomena for comparison with experiment, as this is better done 
with relativistic modeling, but is to provide insight into how the short-range part of the 
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momentum distribution which is important for obtaining saturation of nuclear matter also 
has a role to play in the corresponding scaling phenomena. We shall, however, explore several 
properties of scaling and scaling violations that are observed experimentally, namely, how 
scaling is approached for large momentum transfers q (in the scaling region it is approached 
from above, something most models fail to explain) and whether or not the present model 
is capable of explaining the observed asymmetry found in the scaling function. 

In passing let us draw some comparisons with Deep Inelastic Scattering (DIS) of leptons 
on the proton which has been understood on the basis of the Bjorken scaling law and the 
parton model. There also one observes scaling violations which in the high-energy situation 
are coped with using the so-called evolution equations. Two basic assumptions lie at the 
foundations of Bjorken scaling (see, for example {q]): (1) the highly virtual photon interacts 
with the proton through point-like constituents (the partons) and (2) the partons cannot 
change their momenta during the extremely short time interval available for the DIS process 
and the parton-parton interactions are very weak, a situation referred to as asymptotic 
freedom. When taken into account the latter leads to scaling violations. In contrast, the 
situation of interest in the present work on electron scattering from strongly interacting 
nucleons in nuclei is apparently quite different. The model we use does not display the 
equivalent of asymptotic freedom and yet scaling is quite well obeyed, despite the strength of 
the partonic (nucleonic here) interactions. Of course our model will not be able to account 
for all types of scaling violations, namely those which stem from partonic sub-structures 
(gluons in QCD versus mesons in nuclei via meson exchange currents; the latter have been 
the subject in other studies, for instance, 0, Q). In other words our nucleons (partons) 
are viewed as point-like. The very strong correlations between the nucleons induced by the 
short-range repulsion appear, at least in our model, to lend themselves to a description in 
terms of a mean field framework, in which the nucleons do not interact, apart from Pauli 
correlations. In other words the effect of the hard core is embedded in the modification of 
the nucleon momentum distribution with respect to that of the Fermi gas, still keeping an 
independent-particle model for the system. Thus to the extent that the mean field provides 
a realistic description of nuclei at large momenta, then scaling should occur, as we have 
found. 

This paper, which is closely connected to the research developed in [l2|, outlines the 
model in Sect. [TTl Section UTTl addresses the problems of linking the model to conventional 



3 



perturbation theory and focuses on the Coulomb Sum Rule (CSR), a quantity crucially 
dependent upon the pair correlation function (pcf) or, equivalently, upon the momentum 
distribution n{k) of the nucleons in the nucleus. Indeed the CSR represents one of the best 
testing grounds for the pcf and n{k). In Sect. llVl the fermion propagator, the key element in 
constructing the response of our system to an external probe and hence the scaling function, 
is set up. All the issues connected with scaling and the results we have obtained with our 
approach are collected in Sect. |Vl and then finally in Sect. |Vl] we summarize our main 
conclusions from this study. 



II. THE MODEL 



E AV f dk k^ 1 /■ 



The basic formula we start with reads [13] (we use h = c = 1) 

(27r)3 2m' 

where the factor of 4 accounts for the spin-isospin degeneracy (a summation over these 
variables is of course understood in the second term as well) and A is the particle number. 
Equation ([T]) yields the ground-state energy of the system (actually the energy per particle). 
We apply it to an infinite, homogeneous, non-relativistic ensemble of nucleons, viewed as 
enclosed in a large volume V to be let to go to infinity at the end of the calculation. We 
assume that a two-body force acts between the nucleons and that this is described by the 
potential f (r), where r = \fi — f2\. In Eq. ([1]) C{fi — r2) is the pcf simply related by a 
Fourier transform to the n{k) in our spatially homogeneous system. 

Now, rather than attempting to compute n{k) starting from the potential v{r) adopting 
one of the various many-body techniques available for the purpose (for example a perturba- 
tive one), we assume the momentum distribution to be parametrized as 

nik) = e{kF -k)(l- a-A + e{k - kF)l3ie~^'^^~^\ (2) 

accounting both for the existence of a high- momentum tail in n{k), as suggested by the 
presently available experimental information, as well as standard theory (see, for example, 



a recent treatment of the momentum distribution and spectral function in 



14j). and for 



the Luttinger theorem which guarantees the existence of a Fermi surface for a "normal" 



4 



interacting Fermi system. Of course n{k) (Eq. [2]) should fulfill the constraint 
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^ f 1 - + 3%{Pi + 2/32 + 2) ) = ^h{a, ft, P,) = no, 



(3) 



uq being the system's constant density. 

/^From Eq. (|2]) the pcf is obtained according to the definition 
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and one gets 
C{r) 
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In Eq. (jl]) |\E'o > is the system's ground state and ^7(^1) and ^s{r2) are the fermion fields. 
Furthermore the above formula holds valid for infinite nuclear matter. If we need the pcf 
only for protons, as in the case of the GSR, then the factor ^ should clearly be replaced 
with a factor | in front of the second term on the right-hand side of Eq. ^ and one should 
set riQ = Z/V, Z being the number of protons. As seen in the above equations the pcf 
naturally splits into a direct and an exchange contribution (the first and the second terms 
on the right-hand side of Eq. (jl])). 

In passing we note that by setting a = /3i = /32 = in Eq. (|5]), namely by recovering the 
theta-function momentum distribution of a non-interacting Fermi system, we get back the 
well-known pcf of a Fermi gas 
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which identifies g{r) with . 

Now by inserting Eqs. ([2]) and ([5]) into Eq. ([1]) we obtain the binding energy per particle 
of the nuclear system providing the potential v{r) is known. For the latter we employ a 



schematic model which retains only the basic features of the nucleon-nucleon force since, 
as already emphasized, our aim is not a precise reproduction of the experimental data. 
Accordingly we employ a mixture of a Wigner and a Majorana force, namely 



v{r) = u{r)[l - 7 + 7PJ = < 



+Uo 7 = r <a 

7 = I a < r < 6 , (7) 
6 < r, 



where is the space exchange operator and 7 a parameter varying in the range < 7 < 1. 
In accord with common wisdom we choose for 7 the values indicated in Eq. ([7]) where one 
recognizes the strong short-range repulsion and the moderate intermediate-range attraction 
characterizing the nucleon-nucleon interaction. 

Now all of the elements to compute the behaviour of the E/A versus kp are available. 
We get 

R 4 4-7r 1 /"OO Tin r n 

T = -7^^ / dkk^k) + ^47r / drr' f/o^(a - r) - -^,^(6 - r)e(r - a) (8) 

A uq [Ztt)-^ Zm Jo Z Jo ^ 8 J 

- -^47r- / drr^ Uo9{a - r) + -Vo9{b - r)9{r - a) g^{r) . 
2 4 JO '- 2 -I 

The numerical factors appearing in the potential terms stem from summing over the spin- 
isospin variables of the interacting nucleons; they of course enter differently in the direct and 
exchange contributions to E/A, namely the second and the third terms on the right-hand 
side of Eq. (jH]). 

Before actually displaying the behaviour oi E/A versus kp we have to face the crucial 
problem of fixing the values of the seven parameters (four for the potential v{r), three for 
the momentum distribution n{k)) needed to render our approach predictive. In order to 
tackle this problem we proceed as follows: we start by choosing "reasonable" values. Next 
we compute Eq. (|8]) using these chosen values and repeat the procedure adjusting at each 
step the parameters until they yield E/A ~ -16 MeV and, for the compression modulus, 
K = ~ 14 MeV jisl at the minimum of the curve, that should occur at a value of 

kp which, when inserted into Eq. ([3]), provides the experimental density of nuclear matter, 
namely 0.17 fm~^. 

The three above mentioned constrains (energy, density and compressibility) turn out to 
be fulfilled by choosing 

a = 0.2, /3i = 0.4, 132 = 4: (9) 



and 

f/o = 2.5 GeV, Vq = 53 MeV, a = 0.465 fm, 6 = 2.10 fm . (10) 

With these values the minimum of the curve yielding E/A versus kp occurs at /ci?=1.23 
fm^^, which is obviously different from the value of the pure Fermi gas, namely kp=1.36 
fm~^ . 

The associated E/A versus kp is displayed in Fig. [1] which yields 




= -15.68 MeV, {kF)rain = 1.23 fm-\ {K)min = 13.8 MeV , (11) 

min 



namely the experimental values. 




0.0 0.5 1.0 1.5 2.0 2.5 

kp [fm-^] 

FIG. 1: The binding energy (E/A) versus kp, as given by Eq. 

We also display the momentum distribution n{k) in Fig. [2l Concerning the results of this 
section, it should be clear that the choice of the parameters given in Eq. (fTOj) is far from 
unique, and that these values should be viewed as providing a first orientation on a complex 
problem. Our goal is only to use something representative for the potential in attempting to 
shed light on the connections between saturation of nuclear matter and electron scattering 
scaling phenomena. 

III. THE MANY-BODY CONTENT OF n{k) AND THE CSR 

Next the question arises: which are the Feynman diagrams one has to take into account 
in order to obtain the n{k) given by Eq. ([2])? This is equivalent to asking: what kind of 
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FIG. 2: (color online) Solid (blue online): the momentum distribution n{k) versus k. Also displayed 
(dashed, green online) is the Fermi gas n{k) = 0{kp — k) with kp = 1.36 fm~^. 

correlations among nucleons are responsible for changing the theta-function into our n{k) 
which we assumed to be the true momentum distribution of our system? 

To help in better grasping the relevance of this question it is of importance to realize that 
if we replace in Eq. ([1]) n{k) with 6{kp — k) and C{fi — with Eq. ([6]) we obtain 
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(12) 



namely the result provided by Hartree-Fock (HF) theory ISj, which, as is well known, 
captures the content of the two first-order perturbative diagrams shown in Fig. [3l It is 
instructive to look at the result one gets in HF with the parameters given in Eq. flTOj) for the 
double square well potential of Eq. (j7]). This is displayed in Fig. HI here we see that E/A 
in HF still saturates, although at a far too low density {kp = 1.02 fm~^) and with a far too 
low energy {E/A = -1.56 MeV). 

To prove that by summing all perturbative diagrams one would recover the n{k) of Eq. ([2]) 
(and hence the E/A of Fig. [1]) is clearly an impossibility. In the present study our approach 
is only to explore a simple model while attempting to maintain as hi gh a level of coherence 
as we can. To achieve the latter we first resort to Calogero's theorem [iGj , which states that 
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(a) 



(b) 



FIG. 3: First-order proper self-energy S*-^^; in (a) and (b) are shown the direct and exchange terms, 
respectively. 
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FIG. 4: E/A versus kp in HF theory employing the potential of Eq. ([7]) with the parameters of 

Eq. (uni). 



for a homogeneous infinite system asympotically one should have 



n{k) 



2mnn 



(13) 



where, in the present case. 
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(14) 



In Fig. |5]we compare our n{k) with the Calogero's theorem predictions: we do so in the 
range of momenta starting from the Fermi surface {k = kp = 1.23 fm~^) up to A; ~ 3.6 fm~^, 
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FIG. 5: (color online) Solid line (blue online): the momentum distribution of Eq. ([2]). Dashed line 
(red online): Calogero asymptotic behaviour Eq. (fT3|) . multiplied by a factor 2 x 10""^ (see text for 
the origin of this factor). 

since for larger momenta both distributions become so small that they render the comparison 
meaningless. In Fig. |5]one sees that the two curves are not too different; hence our n{k) and 
v{r), while not exactly the same, nevertheless display an acceptable degree of coherence. It 
should be stated that in carrying out such a comparison, unfortunately, Calogero's theorem 
does not quantify the value of q signalling the onset of the asymptotic regime, and hence 
we have arbitrarily normalized the Calogero's asymptotic momentum distribution in such a 
way to have it coincide with our n{k) at the Fermi surface. 

Concerning the question related to the measurement of n{k), we recall that access to 
information on this is offered by the CSR, which in fact essentially depends "only" upon 
the momentum distribution at least within the context of the present model. Indeed from 
unitarity one has 



where g{r) (see Eq. ([5])) is directly fixed by the Fourier transform of our n{k). 

For sake of completness we display in Fig. |6]the CSR predicted by our n{k). Actually 
this curve was already shown in [l2|, however we revisit it once more here to illustrate how 
the attainement of the asymptotic value (namely one) is postponed by the nucleon-nucleon 
correlations, in particular, the repulsive short-range ones, to larger values of q (= 4.6 fm~^) 




(15) 
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than for the Fermi gas situation (= 2.46 fm ^). 




FIG. 6: (color online) The CSR for a free Fermi gas (dashed line, green online) and for a correlated 
one according to our model (solid line, blue online). 

To iUustrate our second path toward coherence we commence by recahing the ahernative 
model of 17|, which deals with the same issues treated in the present paper, namely the 



sca ling function and momentum distribution, but with a different philosophy. The study 
of [17| is based on the use of PWIA (Plane- Wave Impulse Approximation) and has as its 
starting point the assumption of the factorization of the one-particle exclusive cross section 
(namely the cross section for the process {e,e'N)) according to the expression 

r Hn- 1 PWIA 

. ,,,n',p,,aj =AV-(,..; p. (16) 

In the above if is a kinematical factor, a^^ the eA^ single-nucleon cross section and S{p, £) 
the nucleon spectral function expressed in terms of the so-called missing-energy and missing- 
momentum variables (for their definition see, for instance, isl). We do not dwell here on 
the procedure leading from Eq. (fT6!) to the scaling function: this topic has been addressed 
a number of times in the literature. Suffice it to say that the key step is the integration 
of Eq. ( IT6|) over the variables p and £ in a domain which has been analyzed in the past 
and whose boundaries are set by energy and momentum conservation. The above procedure 
corresponds to passing from the semi-inclusive (e, e'N) cross section to the inclusive (e, e') 
cross section. Indeed one has 



da 
de'dO.' 



a^^(g,a;;p= |l/|,£: = 0)F(g,a;) 



(17) 
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where a^^ is evaluated at the minimum values of £ and p allowed by kinematics. The scaling 
function F{q,u) is easily found to read 

F{q, uj) = 27i [ [ pdpd£S{p, £) (18) 

in PWIA. Now if the integration domain Ti{q,u) in Eq. f fTSj) extends up to encompass the 
limiting value £ — )■ oo then the scaling function can be linked to the momentum distribution 
according to [24 1 



n{p) = / d£S{p, £) . (19) 





Thus, at variance with our method which is based on the coherence between the momen- 
tum distribution and the scaling function, in the sense of having both of them derived within 
the same theory, in [l7] the strategy is to exploit the link between the scaling function and 
the momentum distribution and hence the possibility of extracting information on n{k) from 
(e, e') experimental data. 

IV. SETTING UP THE SINGLE-NUCLEON PROPAGATOR 

Our scheme is based on the central role played by the single-fermion propagator G{k, u) 
for nucleons in the nucleus. In fact, as is well known, knowledge of the latter gives access 



to the E/A according to the expression [18 1 



E = — — hm / dk / c/we*'^" — + uj] TrG(fc, oj) (20) 
2(27r)S->o+i i-oo V2m / ^ ' ^ ^ ^ 

and to the scaling function according to the formula 

F(q, to) = - — -Imn(g, to) . (21) 

m IT 

I — I 

This represents the content of linear response theory and, as in [17], assumes the factorization 
of the single-nucleon cross section. In Eq. ( 12T|) Il{q,u) is referred to as the polarization 
propagator or, in coordinate space, as the density-density correlation function. In field 
theory language it corresponds to a particular choice of the field arguments in the two- 
particle propagator and in momentum space reads 

Il{q,u:) = -'-J-0^^G{k + q)G{k). (22) 
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Knowing G ensures control of n(g, and this, in turn, through Eq. (12T!) . permits the 
determination of the scahng function F{q,u). Also, in a coherent scheme, G{k) should yield 
the correct momentum distribution (namely our n{k)). 

Can such a propagator be derived on the basis of the knowledge of n{k) alone? The answer 
is yes in the simple approach taken in the present work where an infinite, homogeneous 
many-body system of nucleons has been assumed. This kind of mean-field approximation 
possesses a remarkably coherent structure; whether this remains true in a more sophisticated 
many-body framework has yet to be proven. Continuing to work within the context of our 
simple model, as a first step to achieve this goal we rewrite our basic expression in Eq. ([T]), 
exploiting the Faltung theorem of Fourier analysis. For this purpose, we start from 
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with g'^{r) given by Eq. (jl]) and the direct and the exchange potential terms by Eq. 
coordinate space and introducing the Fourier representations according to 
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Then, with the help of Eqs. f l24] - [26|) . we obtain 
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is the single-hole energy displayed in Fig. [71 For comparison, in Fig. [S] the single-particle 
energy 



1 - n{k) (ft) 
n{k) ^'^ 



(29) 
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is shown. Note the discontinuity of ~ 6.5 MeV in both e^.''^ and e^^^ at the Fermi surface, 
the vanishing at large k of e^^'' and the value 



(h) no _ 
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FIG. 7: (color online) Single-particle energy (hole) of our model shown for kp = 1.23 fm ^ which 
corresponds to the saturation density of the system. The dashed line represents the fit given by 
Eq. 



The above single-particle energies yield the poles of the fermion propagator which, as a 
consequence, can then be cast into the form 

n{k) 



G{k,uj) 



^ - + fMO) - ^^(i^ /o" dpp^Mp) /-I dxnilp- k\)] - tv 



(31) 
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FIG. 8: (color online) Single-particle energy (particle) of our model shown for kp = 1.23 fm ^ 
which corresponds to the saturation density of the system. The dashed line represents the fit given 
by Eq. dMI)- 

1 -n{k) 

+ ^ _ [1 _ n{k)][£^ + fMO) - ^^^^^ /o~ dpp^VEip) /-I dxn{\p- k\)] + ' 

This structure of the propagator tells us that in our model the holes exist below, but also 
above, the Fermi surface. Likewise the particles exist above, but also below, the Fermi 
surface. These occurrences clearly reflect the behaviour of our momentum distribution. The 
important point to be stressed, however, is that the propagator in Eq. (pTj) provides the 
correct system energy and n{k\ 

Now, to pave the way to the actual evaluation of 11, it helps to realize that, although 
the expression of the single-particle energies is far from being simple (in fact it cannot be 
expressed analytically), nevertheless its k dependence lends itself to be suitably represented, 
apart from the factor n(/c), by a parabola. This is reminiscent of HF theory. Hence we use 
the following quite faithful representation 

e^''\k) = n{k){A^n) + %)A;^) = n{k)(A^,) + (32) 

with A(^h) = —77.16 MeV and B(^h) = 41.10 MeV fm^, yielding an effective mass m*^^ = 0.50 
m for holes and 

e^^\k) = (1 - n(fc))(A(,) + i?(,)fc2) = (1 - n(fc))(A(,) + (33) 

15 



with = -43.09 MeV and = 22.72 MeV fm^, yielding an effective mass m^^) = 0.91 
m for particles. It is indeed startling to realize how large the impact of our two-body 
interaction in Eg. ([7]) is on the effective hole mass, the generally accepted ratio being actually 



m* /m ^ 0.83 20|. How faithful Eq. fl521) and Eq. fl55]) are in providing the single-particle 



energies, can be garnered from Figs. [7] and [HI 

Finally, the response function of the system is easily derived, 

uj) = -— Imn(g, oj), (34) 

71 

and from the response one immediately obtains the scaling function per proton according to 
q R{q,uj) 

F{q,uj) = —— 35) 

m Z 

= — — At r dkk^n(k) I dx\l - ni\k + q\muJ - e^^\\k + q\) + e^^\k)l 
m uq tx^ Jo J-i 

where the trivial frequency and azimuthal integrations have been performed. The results of 

our numerical calculations of Eq. f l5^ are reported in the next section. Following standard 

practice when discussing scaling we shall show results for the dimensionless scaling function 

f{q,co) = kpX F{q,co) (36) 

which takes on an especially simple form for the relativistic Fermi gas (see, e.g., {2, 3|). 



V. RESULTS 



In Figs. [9] and [TO] the results for the response R{q,co) (Eq. (!M|) ) and scaling function 
f{q, u) (Eq. f l36|) ) are shown versus u for a range of momentum transfers. The response and 
the scaling function for the free Fermi g also displayed for comparison. 

We observe the following: 

• The response and scaling function obtained using our model, as expected, span a range 
of energy loss that extends to larger values than that seen for the Fermi gas model — 
a clear indication of the role of correlations among the nucleons. The widths seen in 
our model are somewhat larger than those of the Fermi gas and the peak heights are 
somewhat lower, both in better accord with experimental data. 

• The peak positions in our model are shifted to higher energy loss than for the Fermi 
gas (see also the discussions to follow). 
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FIG. 9: (color online) The response function of our model (solid line, blue online) and the response 
function of a free Fermi gas (dashed line, green online) plotted versus uj for q = 2kp = 2.46 fm~^ 
up to g = A.bkp = 5.53 fm^^ in steps of 0.5kp. Results for lower values of q peak at lower values 
of UJ. 

• Unlike for the Fermi gas model our R{q,u) and f{q,u) are no longer perfectly sym- 
metric around their maxima. While approximately so, they have tails that extend 
both to higher and lower values of u. However, the degree of asymmetry is not as 
large as what is observed experimentally. 

• Note that while the height of R{q,oj) decreases with q, the height of f{q,u) remains 
constant. 

To investigate the scaling behaviour of our results we follow the usual procedures and 
display f, not versus u as above, but versus the well-known non-relativistic scaling vari- 



Now we indeed see in Fig. [TT] that the scaling functions for different values of q tend to 
group together very closely when displayed versus ipnr, that is they scale. Noting that the 
coalescence in our model occurs at a peak value other than ipnr = 0, it is interesting to 
investigate whether a simple modification of the scaling variable different from the one of 
Eq. ( 1371) can be devised to shift the peak position to zero. One can always do this by 



able [21 1 




(37) 
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FIG. 10: (color online) The scaling function of our model (solid line, blue online) and the scaling 
function of a free Fermi gas (dashed line, green online) plotted versus uj for the same values of q 
used in Fig. EJ 



employing the variable 



, If muj' q 

^nr = Tr[ o 

kp \ q 2 



(38) 



where u' = u — Eshift{q), and where Eshift{q) is a g-dependent energy shift. If one uses the 
simple parametrization 



Eshiftiq) = Eo + E^{q/kp) 



(39) 



with Eq = —17.4 MeV and Ei = 15.9 MeV (see Fig. [T3|) . then almost perfect scaling 
centered about ip'nr = is attained, as seen in Fig. [121 It appears that the impact of the NN 
interactions that have been incorporated in the present model is felt via a g-dependent shift 
in the definition of this new scaling variable, although a direct connection to the underlying 
dynamics is not obvious. As is clear from Fig. [13] similar shifting is seen using relativistic 
mean field theory. One of course should not expect these to be identical, since the present 
model is non-relativistic while the RMF results are obtained using a relativistic model. In 
either case, while the shift is qualitatively what is observed experimentally, it is probably 
somewhat too strong in both cases [22, l23l ]. 

In Fig. [T3I we show the same results as in Fig. [121 but now on a semilog scale. The 
asymmetry, while small, is clearly apparent. More strength is shifted to higher values of 
and, whereas the Fermi gas cuts off abruptly and is only nonzero within the Fermi 
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lAnr 

FIG. 11: (color online) The scaling function of our model (blue online) displayed versus the non- 
relativistic variable iprir- The momentum transfer range is the same as in the previous two figures. 
For reference the Fermi gas result is shown as a dashed curve (green online); clearly, by construction, 
it scales perfectly. 




FIG. 12: The scaling function of our model displayed versus the scaling variable of Eq. ()38p for the 
same values of q used in Fig. [TTl 



cone, the present model produces strength extending to very large and small values of ip'nrj 
in accord with experiment. It is worth remarking that if the high momentum tail in the 
momentum distribution in Eq. ([2]) is set to zero {i.e., (3i is set to zero) then these tails 
extending to large essentially disappear. While setting (3i to zero is not simply setting 
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FIG. 13: (color online) The g-dependent energy shift in Eq. ()39p (solid curve, blue online) together 
with the energy shift obtained in RMF studies of ^^C (dashed curve, red online). 
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FIG. 14: (color online) As for Fig. \T2\ but now on a semilog scale. 

the part of the momentum distribution arising from short-range correlations to zero, since the 
long-range correlations also move some strength from below the Fermi surface to somewhat 
above it, it is very suggestive that in the present model the origin of the tails in the scaling 
function are principally due to the short-range physics, as is often assumed to be the case. 
Interestingly, the position of the peak is largely unaffected by "turning off' the high-fc part 
of the momentum distribution, suggesting that the peak position is not strongly correlated 
with the short-range physics, but rather with the long-range physics. 
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FIG. 15: The scaling function versus q for [V'nrlo = —0.6, -0.4 and -0.2, i.e., in the scaling region. 

In concluding this section we address the problem of how the scaling regime is approached. 
We do so using the scaling variable [V'nrlo) where, to make this closer to what has been used 
in analyses of experimental data [1] a constant energy shift Eghift = 30 MeV has been chosen 
in using Eq. fl38l) : this is indicated by the subscript "0". We display / as a function of q 
for three values of the scaling variable in the scaling region, i.e., to the left of the QE peak. 
The curves are shown normalized at g = 2kF, namely, the ratio 



is displayed. From Fig. [15] it clearly appears that the scaling regime is approached from 
above, an occurrence which is qualitatively in accord with experimental findings. 

VI. CONCLUSIONS 

In the present study we have developed a model centered around an assumed form for the 
momentum distribution of nucleons in nuclei. The momentum distribution has been chosen 
to reflect current understanding of how dynamical effects underlying the nuclear many- 
body problem lead to a form for n{k) with both low-k components coming from long-range 
interactions and an extended tail at high-A; arising from short-range interactions. We have 
restricted the scope of the study to infinite, homogeneous nuclear matter, have employed 
only point nucleons, have assumed a strictly non-relativistic model (although in future work 



piq) = fiq, K.]o)//(2/.F, KM 



(40) 
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we hope to extend the scope to relativistic modehng), and have restricted our attention to 
the longitudinal electromagnetic response, for the present. 

Working with this as a basis we have developed the formalism in two different directions, 
maintaining as much consistency as possible. First, we have devised a single-particle Green 
function that leads to the know saturation properties of nuclear matter and to realistic 
particle and hole single-particle energy spectra. Second, we have taken the same Green 
function to obtain the density-density polarization propagator and, through its imaginary 
part, have obtained the longitudinal electron scattering response function R{q, ou) and the 
scahng function f{q,io). For the latter we have explored several aspects of scaling and of 
scaling violations. 

We find that scaling is quite well respected, despite the strength of the NN interactions 
implicit in the problem. There are seen to be some scaling violations, for instance, those 
observed as shifts of the quasielastic peak positions as functions oi q — indeed these are 
also seen in relativistic mean field theory. The shape of the scahng function is observed to 
be more spread out than is the Fermi gas scaling function and the former yields a somewhat 
lower peak height than the latter, both in rough accord with experiment. Where it is large 
the scaling function in our model is somewhat asymmetric; however, it is not enough so 
to agree with experiment. When the tails of the scaling function are examined in detail 
the asymmetry is more apparent and, indeed, the strength at both very large and very 
small values of the scaling variable is significant, in accord with experiment. If the high- 
k contributions to the momentum distribution are "turned off' then these tails disappear, 
suggesting that their origin lies in the part of the momentum distribution arising from short- 
range correlations. The position of the peak, however, appears to be due to long-range 
physics. 

In summary, the present study has demonstrated that a high level of consistency can be 
maintained in simultaneously representing both the saturation properties of nuclear matter 
and the scaling properties of the longitudinal electron scattering response. 
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